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Abstract Critical phenomena of ferromagnetic tran- 
sition in colossal magnetoresistance manganites are the- 
oretically studied. Concerning the critical exponents for 
this transition, there still remains controversy among ex- 
perimental results. In order to clarify intrinsic physics 
of the manganites through a comparison with theoreti- 
cal prediction, we investigate the critical phenomena of 
double-exchange models by using finite-size scaling anal- 
ysis on unbiased nunrerical results. As a result, we show 
that the critical exponents of the ferromagnetic transi- 
tion of the three-dimensional double-exchange model is 
consistent with those of the Heisenberg model, but are 
distinct from the mean-field one. 



1 Introduction 

Ferromagnetic transition in doped manganese oxides has 
revived interests since rediscovery of colossal magnetore- 
sistance phenomena, jy,^ Recent refinement of exper- 
imental technique and improvement of sample quality 
make it possible to discuss critical phenomena of this 
transition in detail. Particularly, for the purpose of cap- 
turing intrinsic physics in these materials and promoting 
a better understanding of the colossal magnetoresistance 
phenomena, the critical exponents and the universal- 
ity class have been investigated by various experimen- 
tal methods. However, the estimates of the exponents 
are scattered and have not converged yet. For instance, 
the neutron scattering gives Heisenberg-like universality 
class, HE! an d some d.c. and r.f. magnetization mea- 
surements give mean- field-like one. (Spl It is strongly 
desired to solve this controversy. 

The purpose of this work is to shed light on this prob- 
lem from a theoretical viewpoint. We study the univer- 
sality class of the ferromagnetic transition due to double- 
exchange (DE) mechanism, which is believed to be a ma- 
jor driving force to cause the metallic ferromagnetism in 
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doped manganites. We compare the results with experi- 
mental data, and discuss the origin of the experimental 
controversy. 

The DE mechanism itself was proposed more than 
half-century before, W and has been studied intensively. 
Nevertheless, the thermodynamics has not been fully un- 
derstood thus far. Large spin fluctuations and the strong 
interplay between charge and spin degrees of freedom 
make it difficult to describe finite-temperature proper- 
ties quantitatively, especially near the critical point of 
the ferromagnetic transition. In order to give a quan- 
titative description of the critical phenomena, it is in- 
dispensable to take account of the large fluctuations. 
Approximational methods like the mean-field approxi- 
mation are known to fail. 0] 

As one of promising candidates to study this prob- 
lem, in this paper, we employ a numerical calculation 
which includes the fluctuations properly and is free from 
uncontrolled approximations. Recent development of al- 
gorithm and computational environment enables us to 
perform large-scale unbiased calculations and to describe 
the finite-temperature properties in the thermodynamic 
limit. We analyze the numerical results which are ob- 
tained by a newly developed Monte Carlo (MC) method 
in the next section 0. In Sec. ||, we discuss the results 
and compare them with experimental results. Section |J 
is devoted to summary. 



2 Numerical Results 

We develop and employ a new MC algorithm, which is 
called the moment-expansion MC method. ||] By using 
the moment expansion of the density of states which 
is efficiently performed on parallel computers, this new 
method reduces the computational time remarkably com- 
pared to the conventional MC technique. 

We study DE models with a single conduction band 
(which stands for Mn e g bands) which couples to local- 
ized spins (in Mn ti g localized levels). The Hamiltonian 
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Fig. 1 Temperature dependence of the magnetization in the 
thermodynamic limit. The curve in the figure is the least- 
squares-fit toma (T c — T)' 3 with the Heisenberg exponent 
P — 0.365. The inset shows the system-size extrapolation of 
the spin structure factor. 



is given by 
H 






(1) 



where the first term describes the electron transfer be- 
tween nearest-neighbor sites, and the second one is for 
the Hund's-rule coupling ( J H > 0) between the conduc- 
tion electron spin <x and the localized spin S. For sim- 
plicity, hereafter, we consider the strong coupling limit, 
Jh — > oo, and the classical localized spin, |S| — > oo. The 
following calculations are done at quarter-filling. 

First, we study a three-dimensional DE model with 
Heisenberg spin symmetry, which is realistic for com- 
parison with experiments. M We analyze the MC data 
by assuming the universality class of Heisenberg models 
with short-range interactions. Figure [l] shows the tem- 
perature dependence of the magnetization in the ther- 
modynamic limit. The data are obtained by system-size 
extrapolation of the spin structure factor as shown in the 
inset. The magnetization data scale well with the fitting 
m ex (T c — T)@ with the Heisenberg exponent f3 — 0.365. 
Figure |2| shows the finite-size scaling for the spin struc- 
ture factor of the ferromagnetic component 5(0). The 
scaling relation is given by 



5(0)1 
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where L is the linear dimension of the finite-size clusters, 
/ is the universal scaling function and e = (T — T c )/T c . 
In Fig. |2|, we assume the Heisenberg exponents, v — 
0.705 and r\ = 0.034. The MC data for different L and T 
appear to follow a universal function. We have tried the 
mean-field exponents f3 = v = 0.5 and r\ = also, and 
we confirm that the Heisenberg exponents give better 
scaling fits than the mean-field ones. 
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Fig. 2 Finite-size scaling plot for the spin structure fac- 
tor in the three-dimensional model. Critical exponents in the 
Heisenberg universality class are assumed. 
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Fig. 3 Finite-size scaling plot of the spin structure factor in 
the two-dimensional model. Critical exponents of Ising uni- 
versality class are assumed. 



Next, in order to perform more intensive and pre- 
cise MC samplings in large-size systems, we study a 
simplified model of eq. (|l|), that is, a two-dimensional 
DE model with Ising spin symmetry. [|10f The localized 
spins are Ising variables which take S = ±1. The model 
is much simplified, however, it still include the essence 
of the DE mechanism, that is, the kinetics of electrons 
is strongly correlated with spin degrees of freedom and 
energetically favors ferromagnetism. We analyze the MC 
data on spin structure factor by finite-size scaling analy- 
sis as in Fig. 0. Figure shows the result on the assump- 
tion of the exponents for the Ising spin system in two 
dimensions, v = 1 and 77 = 0.25. All the MC data are 
scaled well with a universal function. In this case also, 
we confirm that the mean-field exponents, v — 0.5 and 
?/ = 0, do not give a satisfactory scaling. 

This obvious result in the two-dimensional model 
strongly supports the previous conclusion in three di- 
mensions, i.e., the critical exponents in the DE model are 
consistent with the universality class of the Heisenberg 
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models with short-range interactions. The estimates of 
the exponents are distinct from the mean-field ones. These 
results are discussed and compared with experimental 
results in the next section. 



3 Discussion 

The DE interaction has a distinguishable property com- 
pared to ordinary exchange interactions in spin systems. 
Effective ferromagnetic interaction is induced by the ki- 
netics of electrons which favor extended states with fer- 
romagnetic spin background to gain the kinetic energy. 
If we integrate out the electron degrees of freedom to de- 
scribe the action as a function of spin configurations, it is 
necessary to introduce effective long-range two-spin in- 
teractions as well as multiple-spin interactions which de- 
pend on sizes and shapes of ferromagnetic domain struc- 
ture. The range of the interaction is determined to mini- 
mize the total free energy for charge and spin degrees of 
freedom. As the system approaches the critical point, the 
magnetic domain structure fluctuates strongly. Thus, it 
is highly nontrivial how the DE interaction is renormal- 
ized, whether it is renormalized to a short-range one, or 
the long-range and the multiple-spin interactions become 
relevant to cause the mean- field- like transition through 
suppression of fluctuations. 

The results in the previous section indicate that the 
universality class of the ferromagnetic transition in the 
DE model is consistent with that of models with short- 
range interactions. It is distinct from the mean-field uni- 
versality class. Thus, we conclude that the DE interac- 
tion does not become effectively long-range, but is renor- 
malized to be short-range near the critical point. 

As mentioned in Sec. |], experimental estimates of 
the critical exponents are widely scattered including the 
Heisenberg-like (short-range) and the mean-field-like (long 
range) ones. Our numerical results indicate that the short- 
range Heisenberg universality class should be observed 
in real materials when the ferromagnetic transition is 
ascribed to the DE mechanism as a major origin. 

In order to solve the experimental controversy, the 
following points should be considered carefully. The first 
point is sample quality. For precise estimates of the criti- 
cal exponents, it is crucial to determine the critical tem- 
perature T c precisely. Sharpness of the transition as well 
as the value of T c depends strongly on homogeneity of 
samples. The second point is the range of the critical 
region. The exponents may depend on the region where 
they are estimated, particularly when different interac- 
tions compete with each other. Generally, it is difficult 
to determine true critical region, but the dependence of 
the exponents on the critical range should be tested. The 
third point is that, in real materials, there may be other 
elements which are neglected in the simple DE model (Q), 
for instance, orbital degrees of freedom, Coulomb repul- 
sion between conduction electrons, and electron-phonon 
coupling. 



We comment on the final point. We believe that these 
additional elements neglected in model (M do not alter 
the short-range universality class of the ferromagnetic 
transition because they will not suppress the fluctua- 
tions. The spin fluctuations may become more strong 
through these additional channels. At the same time, 
the electrons become more incoherent by these addi- 
tional elements. This makes the effective magnetic inter- 
action more short-range since the incoherence shortens 
the range of the effective electron transfer. These con- 
siderations support our conclusion that the short-range 
Heisenberg universality class will be found in experi- 
ments if the transition is purely magnetic. 



4 Summary 

We have investigated the critical phenomena of the fer- 
romagnetic transition by the double-exchange mecha- 
nism in order to discuss the experimental controversy 
on the estimates of the critical exponents. We have esti- 
mated the critical exponents by using the finite-size scal- 
ing analysis on the Monte Carlo data. The obtained ex- 
ponents suggest that the universality class of the three- 
dimensional double-exchange model is the same as that 
of the Heisenberg models with short-range interactions. 
The double-exchange interaction is renormalized to be 
short-range when the system approaches the critical point. 
This theoretical result provides a useful reference to con- 
sider the experimental controversy. We have discussed 
possible origins of the discrepancy among experiments. 
Although our conclusion is favorably compared with the 
neutron scattering data, further experimental studies are 
desired to settle this controversy. A conclusive determi- 
nation of the critical exponents and the universality class 
will help us to understand the intrinsic physics of doped 
manganites and of the colossal magnetoresistance phe- 
nomena. 
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